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Abstract. - We develop a non-equilibrium Ginzburg-Landau type theory of the far from the 
equilibrium dynamics of superconductors in a fluctuational regime and apply our approach to 
quantitative description of a superconductor island in a stationary nonequilibrium state. We 
derive the effective temperature of the nonequilibrium state and find fluctuational contributions 
to the magnetic susceptibility showing that it becomes a singular function of \/V — Vc, where V 
is the external drive and Vc is its "critical" value at which the nonequilibrium phase transition 
takes place. 
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Time-dependent Ginzburg-Landau equation (TDGL) 
Successfully describes weakly non-equilibrium dynamics of 
the systems experiencing a second order phase transition, 
ificluding the ferromagnetic, superconducting, and the su- 
perfluid transitions, to name a few [THlOj. in the vicinity 
qf the critical point. A marked progress has been recently 
achieved [5] in formulation Keldysh technique-based ap- 
proach aimed at extension of TDGL onto strongly non- 
equilibrium situation, see, e.g. the review Ref. [11]. Yet 
constructing a theory of far from the equilibrium dynam- 
ics for the GL systems remains a major challenge of the 
ijonequilibrium statistical physics. 

' A nonequilibrium extension of the Ginsburg-Landau 
theory (NGL) requires, in principle, the details of the un- 
derlying microscopic physics. The latter enter the the- 
ory through the quasiparticle density of states and relax- 
ation rates that appear as parameters in corresponding 
kinetic equations (KE) for quasiparticle distribution func- 
tions which, in their turn, depend on the order parameter 
A [3115 . However, as we show below, one can construct 
a phenomenological non-equilibrium theory in a critical 
region without invoking the details of the microscopic be- 
havior of a low-symmetry phase, making use of the sym- 
metry considerations in a spirit of Ref. [T2| . 

In this Letter we develop a theory of the far from- 
equilibrium fluctuation effects in superconductors gener- 
alizing a standard description of superconducting fluctu- 
ations [3 US] to the far from the equilibrium state and 



derive a closed GL-like equations set, describing quan- 
titatively far from the equilibrium dynamics of fluctu- 
ations. In particular, we calculate the nonequilibrium 
fluctuation contribution to magnetic susceptibility and 
find an analytical expression for the effective temperature, 
Tcff = T cosh^{V/4T), depending on the bias, V. As an 
illustration of the proposed general technique, we discuss 
the dynamics of a superconducting granule connected to 
reservoir via disordered normal wires. 

In the framework of a phenomenological theory of a 
second order phase transition, the properties of the low- 
temperature phase near the transition are determined by 
the free energy functional written as an expansion in the 
order parameter, see, e.g, Refs. [12l[T4|. In an equilib- 
rium, the density matrix space M is parameterized by 
the temperature and the transition occurs at T ~ T^, the 
corresponding dimensionless parameter of the Landau ex- 



pansion is |Tgl 



the units where h 



|r - Tcl/Tc < 1 (we win use 
= e = 1 throughout the paper), 



where is the Ginzburg-Landau relaxation time in an 
equilibrium. In a general nonequilibrium case the transi- 
tion extends over some surface in M (equilibrium density 
matrices form a zero measure subspace of M) . 

We consider a system where the excitations and their 
kinetics are well defined. Then a nonequilibrium GL the- 
ory has the form similar to that of the equilibrium one, 
provided the distribution functions of the excitations are 
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stationary. Following the general recipe [TT] for treat- 
ing an out of the equilibrium system, one is to use the 
Keldysh partition function, instead of the conventional 
partition function for the equilibrium case, and, accord- 
ingly, the Keldysh action for the order parameter replaces 
the GL free energy in equilibrium. The parameters of the 
noncquilibrium GL-expansion are functionals of the exci- 
tation distribution functions. The closeness of the system 
to the (non-equilibrium) phase transition surface is deter- 
mined by the dimcnsionless parameter |TGLrcff|~^ ^ 1, 
where Tofi is the noncquilibrium energy scale replacing T 
in an out of the equilibrium state and Tql is a general, 
noncquilibrium, GL relaxation time replacing t^^^ . This 
inequality is the necessary condition for applicability of 
our approach. Note that the existence of Toff docs not 
imply the local equilibrium form of the excitation distri- 
butions. 

The coefficients of the noncquilibrium GL functional be- 
have differently as compared to those of the equilibrium. 
In particular, the coefficient at the forth order term in the 
order parameter. A, can even change the sign at large driv- 
ing forces. This signals the onset of an instability of the 
NGL-equations solution and means that the NGL func- 
tional should be expanded to higher orders with respect 
to A. In the context of superconductivity, the equilibrium 
GL-expansion has the usual form jl4] : 
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Fig. 1: (color online) a) Density matrix space, b) An exem- 
plary system: superconducting granule weakly connected to 
the reservoirs, c) The GL nonequilibrium relaxation time as a 
function of the applied voltage V at the resevoir temperature 
r = 0.7Tc. [The units of t^I are chosen to match ln(r/rc) 
in an equilibrium]. The lower inset illustrates that T^s differs 
from V . The upper inset shows biV) and illustrates the differ- 
ence between the nonequilibrium and equilibrium behaviors; in 
the latter case h is the temperature independent constant. 



where L'^ = id^ + (r^.^V - Ddl Ul^Y' = (T - 
Tc)a, Za = J DAD A* c^p{-n[A, A*]} is the partition 
function, and u is the density of states (DoS) at the Fermi 
shell. In a disordered metal model a = tt/ST, a = S/tt, D 
is the diffusion coefficient and b = 7C,{3)/tt^ [T]. 

The order parameter in the Keldysh space has two com- 
ponents, corresponding to the lower and upper brunches 
of the Keldysh time-contour. To simplify the structure of 
the action, it is convenient to use the rotated basis and in- 
troduce "classical" (Ai) and "quantum" (A2) components 
of the order parameter [half-sum and half-difference of the 
order parameter values at the lower and u|)per brunches 
of the Keldysh time-contour] [9lfT5) . Thus A = (Ai,A2)^. 
The Keldysh analog of the partition function Za is 



Za = J D'lD'l*exp{iSA[^,^*]}■ 



m 



The average order parameter should be identified with 
(Al), while (A2) = 0, where averaging is based on Za- 



= J D'^D'^*[...]exp{^SA[^,^*]}■ 



i3) 



Here we used the fact that in the absence of the quantum 
components of external source-fields, Za = 1, this is a 
manifestation of the causality principle. The saddle point 



equation 



SSaI^,^* 
SA* 



= 0, 



(4) 



at the manifold A2 ~ A2 = produces the nonequlibrium 
generalization of the GL equations. 

— ^ — ^ 

Near the phase transition we can expand iS a[A, A*] 
over A. This expansion over quantum components of the 
order parameter should be performed up to the first order, 
see Eq.(|3]). Carrying out the microscopic calculation based 
on the Keldysh functional representation of the BCS the- 
ory fl6| in a form of the nonlinear cr- model we find that the 
shape of S'AiA, A*] above and below [the restrictions are 
discussed] the phase transition can be written in a form 
similar to that of ([T]): 



5a ~ Tr 



L- 



Jeff 



(5) 



where Tr means the trace with respect to times and the 
integration over coordinates, while dx is the Pauli matrix 
in the Keldysh space. We show below that Eq.© can be 
derived as well from the symmetry considerations similar 
to those of Ref. P]. 

The microscopic Keldysh nonlinear a-model calculation 
expresses GL coefficients a and b through the gauge invari- 
ant electron and hole distribution functions [3] , F^. and 
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[or = {Fe + fh)/2 and Ft = {F^ - ^h)/2], where in an 
equilibrium, -Fe(h) = tanh(e/2T): 



\imd,FL{e) = 

e-i-O 8 J off 



4 e-i-0 



tY{t)dt, 



(6) 
(7) 



where 27rf(t) = de cxp[~iEt]Fi^(^E) / {E + iO). One 
can see straightforwardly that in the equilibrium limit 
Eqs.®-© reproduce conventional values a = n/ST and 
b = 7({3)/n^ respectively. Formulas © and ([7|) for the 
coefficients of the far from equilibrium GL functional con- 
stitute the main result of our work. 

A structure of L in a Keldysh space is determined by 
the causality principle [TT] and by the comparison with 
the standard form of GL expression given by Eq.(IT]) 







L 



(8) 



where R{A) and ii'-subscripts denote the retarded (ad- 
vanced) and Keldysh propagators respectively. The form 
of the i?(A)-component can be found from the correspon- 
dence between the imaginary and real time representations 
that can be transformed into each other by the Euclidean 
rotation: 



(9) 



where {^p,A) is the electromagnetic potential [T7] ]. 

The Keldysh component of the fluctuation propagator 
is expressed via the additional temperature scale T*: 



[L 



(-1)1 



AiT* 



(10) 



where T* ~ [1 — Fc(0)Fh(0)]Tcff is determined from the mi- 
croscopic calculation and is determined by (provided 
the time reversal symmetry is not broken): 



(11) 



where i^(^'(e) = a/2e. Here Fj^\e) is the distribution 
function at the phase transition surface. 

The current can be found analogously to Ref. [5] by 
adding to the action the quantum source in a form of the 
vector potential, A2, in addition to the (classical) external 
fields potentials {f, A) discussed above, and then varying 
the term proportional to A2 in Sa- Thus, yields the su- 
percurrent as, 

j, = TTei^D\Ai\^draig{Ai)/2T,tf. (12) 

The conditions that electromagnetic fields and quasi- 
particle distribution functions are stationary and weakly 
depend on coordinates on the scale of the Cooper pair 
size, ^ = VDtgl, ensure the applicability of the NGL- 
functional. We verify that our highly NGL equations re- 
duce to the standard weafcZy-noncquilibrium form in the 



case where the deviations from the equilibrium are small. 
Indeed, if F differs slightly from tanh(e/2T), one has to 
use the equilibrium values for a and b and keep nonequilib- 
rium Tgl [since the nonequilibrium corrections are of the 
next order if the deviation from the equilibrium is small] . 
This means that our theory is restricted to the neighbor- 
hood of the intersection of the phase transition surface 
with the equilibrium density matrix subspace, see FiglTJ 
Rewriting Eq.dni) as l/iT^^T^s) = a{T-T,)+a /jFL(e)- 
tanh(e/2r)]/2e wc recover the weakly- noneciuilibrium ver- 
sion of the GL theory . 

Most of the quantities related to superconducting fluc- 
tuations are the singular functions of Tq^^^ in the equi- 
librium [13]. When we move out of the equilibrium, 
the fluctuation corrections are parameterized by the 
nonequilibrium Tgl and depend on Toff- Recently zero- 
dimensional superconducting fluctuations and fluctuating 
diamagnctism in the lead nanoparticles were experimen- 
tally investigated, see, e.g., Ref. [19]. Motivated by the 
experiments we flnd, as an example, the nonequilibrium, 
V > Vc, fluctuation contribution to the magnetic moment 
of a small (of the size L ^) superconductor: 



AI = -2T^sTc,i.HDL^ 



(13) 



where H is the magnetic field and 77 = 1/10 for a spheri- 
cal island. Nonequilibrium fluctuation corrections to other 
quantities, c.g, related to the diffusion propagator contri- 
butions to Z [and the Langcvin noise corrections with the 
correlator proportional T*], will be presented in |18| . 

Now we sketch a general procedure for calculating 
fluctuation-related quantities. In an out of the equilib- 
rium state, one uses the Keldysh real time partition func- 
tion instead of the conventional thermodynamic partition 
function [3]: 



Z= D[A,'i>,'i>A,A*]exp{iS}, 



(14) 



where the Grassman flelds, 5'. 5*, describe the 
fermion (superconductivity-related) degrees of freedom 
(on Keldysh contour), A = (<^, A) and S'[^, 5', §A, A*] 
is the microscopic Keldysh action of the system. The 
Ginsburg-Landau expansion of the effective thermody- 
namical potential is an example of the so-called low energy 
field theory, i.e. a theory, where the order parameter fields 
change negligibly on the microscopic scales, e.g., the lat- 
tice constant. The low energy theory in Keldysh formalism 
appears after integrating out the high energy part of the 
fields. The resulting low energy effective action consists 
of the three parts. The first one, Sa, describes quantum 
dynamics of the A-field, the second part generates the ki- 
netic equations for the excitations. It looks schematically 
like Tt[z o (KE)], where the dynamic variable z is closely 
related to the anti-Keldysh component of the Q-matrix in 
the nonlinear cr-model formalism. Integrating over z we 
obtain the functional i5-function ensuring that the distri- 
bution function obeys the kinetic equations. The third 
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part Sa describes electromagnetic fields. The variation of 
the effective action over Sa produces Maxwell equations. 

The fluctuations of the A-field enter the collision inte- 
grals of the kinetic equations and the collisionless terms 
[the fluctuation renormalizations of the KE-cocfficicnts] , 
while (Al) enter the (nonlinear) kinetic equations as ex- 
ternal fields |4ll8j . The A- fluctuations in KE contribute to 
the fluctuation corrections to the kinetic coefficients [T^. 
The (Ai)-tcrms in KE are important while differs es- 
sentially from tanh(e/2T) only at small energies, e ^ (Ai). 
But if quasiparticles are excited in the wide energy range 
above the gap, (Ai) < |e| < Tcff, then (Ai)-terms in KE 
induce small, ~ o((Ai)/rcff) ^ 1, the perturbation of 
and the subleading (Ai)/T'cff ^ 1] terms in com- 
pared to the terms given in Eq.©- 

The phenomenological Landau theory predicts (tql)~^ 
to depend linearly on \T — Tc\ and Itg^ITc ^1- In the 
nonequlibrium state the role of {T ~ Tc) is taken by some 
functional of the electron and hole distribution functions, 
which characterizes the effective "distance" from the phase 
transition. We expect that t^^ is a linear functional of 

SFe{h) = -^L(T) ~ ^L(T)' Otherwise the contribution pro- 



L(T) 



portional to (T — Tc)"^ to r^^ would have appeared in an 
equilibrium. We write thus 



(15) 



where the kernels K are some functions of the energy. In 
an equilibrium Ft = and .Fl — tanh(e/2T). Then we re- 
produce the equilibrium value of t^^^ = a{T — Tc) choosing 
K^^\e) = a/2e. We consider the system invariant under 
the time reversal symmetry. So we should choose K^^^ = 
because otherwise this term would give the unnatural con- 
tribution to changing its sign when, e.g., we reverse 
the direction of all currents in the system. 

The important question is how this formalism describes 
the phase transition interface in the density-matrix space, 
see FiglT] The parameter r^T^ should not depend upon 



the choice of Fj^^ , thus F-^"' and F^' ' belonging to the 
interface should satisfy the relation 



AO) 



(0' 



de 
Ye 



Ff\e)~Fr{e) 



io)i 



0, 



(16) 



where Wd is the Debye energy. 

The microscopic derivation of the GL action shows that 
Eq. can be interpreted as the integral representation 
of the electron-phonon interaction constant. A, giving the 
BCS-superconductivity: 



(0) 



^tanh(6/2Tc) = -i^, (17) 



with C = 0.577. 



where Tc = 2^u}o/Tre~^/'^^ and 7 =: 
being the Euler constant. 

In the Fourier space -^^(^a) ^ ■^^'^ + '''gl + Dq^. The 
Keldysh component ^' in an equilibrium should sat- 
isfy the relation following from the fluctuation-dissipation 



theorem (FDT) [3]: 



(18) 



where Bi^ = cothw/2T is the equilibrium distribution 
function of the complex A-field. The similar relation holds 
for the out of the equilibrium state where the gradients of 
B with respect the "the center of mass" Wigncr trans- 
formation variables arc irrelevant, which is the case we 
consider. The main (infrared) frequency scale of the Lan- 
dau theory is r^T^ and Dq^ ^ ~ ■ In an equilibrium 
Tgi^ (T-Tc) < Tc and, therefore, 2iu}B^ AiT^. In the 
out of the equilibrium state we should choose smaller 
than any relevant energy scale of B^. Then we can also 
replace 2a; i?^ by 



lim 2iuB^ = T* 



(19) 



and consider T* as the second effective temperature. So 
[L("1)]k = 4ir*. (20) 

A recipe for constructing the stationary nonequilibrium 
distribution function is sketched in FigUlD. The supercon- 
ductor is connected to the electron reservoirs through the 
wires with the normal resistances R and R' [e.g., the disor- 
dered quasi-lD normal metal wires]. The reservoir biases 
and temperatures arc, in general, different. If the wire 
resistance satisfies the relation, R 3> l/La,^ [L and ctn 
are the S-diameter and the normal conductivity of the S- 
material] then the current coming from the wires spreads 
over the island that ensures the weakness of fc(h) (f) gra- 
dients. The supercurrent js is much smaller than the crit- 
ical current if ^reff/(i(Ai)) La^R [the normal current 
- supercurrent conversion at the NS interface is consid- 
ered in Ref. fl8| using the generalized boundary condi- 
tions [20j . which take into account the proximity effect]. 
We neglect the Coulomb blockade effects considering the 
transmission probability (per channel) between the island 
and the wires to be close to unity. Normal wire - su- 
perconductor island structures used in the experiments 
PT| satisfy most of these conditions. Solving the nor- 
mal state KE, Dd^Fc(h) = 0, we find that on the island, 

^e(h) = E„Pntanh[(i; - Vn)/2T„], EnPn = 1, whcrC Vn 

is the voltage at the terminal n and Pn is determined by 
the resistivity of the wires [the exemplary system size is 
smaller than the quasiparticle inelastic length scale]. For 
the two terminal case, p2 = R/Ry,, where i?5] ~ R + R' . 
Similar structure of p„ holds for a multiterminal case. 

To illustrate the developed approach, we consider an 
exemplary system with P2 = Pi = 1/2. Then 



F 



0(h) 



tanh 



E±V/2 
2T 



tanh 



EtV/2 



2T 



(21) 



p-4 



Nonequilibrium Ginzburg-Landau-type theory 



where Sf = 0^,T*^Tcff, and 

Teff = Tcosh2(l//4T), (22) 
r-J = aTeff {RcH_i+„ + 21n2 + hr(r/Te)} , (23) 

6 = -;^RcM'(2)(l/2 + zu). (24) 

Here u = V/inT, H is the Harmonic number and is 
the Digamma function. Importantly, b is very sensitive to 
the degree of nonequihbrium, see FigU] it changes sign at 
u w 0.3 remaining negative at larger u that signals of the 
possible instability [in agreement with [lOj ] and requires 
keeping ~ A^-tcrms in the action [TS]. While u < I, r^T^ ~ 
[Ou^ + ln (T/Tc)] Q;Tcff with = 7C [3] ^ 8.4. In FigUl 
(A) - VVc - V, with Vc ~ 47rT^ln(Tc/T)/6i. Taking, 
e.g., T = l.STc ill the reservoirs and V ~ 2Tc we get: 
Teff/T - 0.87, tgl/t^''^'' ~ 0.69 and then the fluctuation 
susceptibilities, see Ea. (fT3| : x/x^*"*^ ~ 0.6. 

To conclude, we have constructed the nonequilibrium 
GL theory on the symmetry grounds under the condition 
that the kinetics of the high-symmetric phase is estab- 
lished. The coefficients of the nonequilibrium GL func- 
tional, which are the constants in an equilibrium, become 
strongly dependent on the external drive in an nonequi- 
librium state. In particular, the coefficient at the forth 
order (in the order parameter) term can change its sign 
at large driving forces; this would signal the onset of the 
instability which requires the higher order expansion. The 
energy parameter Tes replacing the equilibrium tempera- 
ture T, is now a nonlinear function of the bias voltage. 
We have demonstrated that the fluctuation corrections to 
observable quantities, e.g., to the magnetic susceptibility, 
in a superconducting island get strongly renormalized and 
become the singular functions of \/V — Vc when out of the 
equilibrium rather than being functions of \/T — Tc in an 
equilibrium. Accordingly, the order parameter vanishes 
like \/Vc — V in the out of the equilibrium state replacing 
its \/Tc — T-dependence of an equilibrium state. 
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